Metastable quantum phase transitions in a periodic one-dimensional Bose gas: 

II. Many-body theory 
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Wc show that quantum solitons in the Lieb-Liniger Hamiltonian are precisely the yrast states. We 
identify such solutions with Lieb's type II excitations from weak to strong interactions, clarifying 
a long-standing question of the physical meaning of this excitation branch. We demonstrate that 
the metastable quantum phase transition previously found in mean-field analysis of the weakly 
interacting Lieb-Liniger Hamiltonian [Phys. Rev. A 79, 063616 (2009)] extends into the medium- 
to strongly interacting regime of a periodic one-dimensional Bose gas. Our methods are exact 
diagonalization, finite-size Bethe ansatz, and the boson-fermion mapping in the Tonks-Girardeau 
limit. 
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I. INTRODUCTION 

Exactly solvable quantum systems P, Q are now within 
reach of experiments. This is best accomplished in highly 
controllable systems, such as ultracold quantum gases Q, 
because one has precise control over the effective dimen- 
sionality, so that one and two dimensions can be studied 
for a wide range of interactions, both repulsive and at- 
tractive, from the weakly to the strongly interacting, over 
seven orders of magnitude Moreover, these systems 
are well insulated, where there is only negligible exchange 
of energy and particles with the environment, and thus 
suitable for the study of the metastable quantum phase 
transitions in excited states as well as ground-state 
quantum phases. 

In this article we investigate the many-body ground 
and excited eigenstates of a periodic one-dimensional 
(ID) Bose gas [6[ under an external rotating drive going 
beyond the mean-field regime. Such a geometry has been 
realized in experiments |7Hlll] from the weakly interact- 
ing condensate regime to the strongly interacting Tonks- 
Girardeau (TG) limits [l^ - [T^ . In our previous analyses 
we showed that the average angular momentum of weakly 
repulsive bosons in a one-dimensional ring undergoes a 
quantum phase transition (QPT) in the metastable states 
as a function of interaction and rotation [l5j . In the 
mean-field theory this phenomenon is intuitively under- 
stood in terms of bifurcation of stationary excited-state 
energy branches of the plane-wave state propagating on 
the ring, and of localized soliton trains |16[ ■ Each excited 
state has a denumcrably infinite number of bifurcations 
from the plane wave to a state containing one or more 
gray or dark solitons; each such bifurcation corresponds 
to a QPT. Formally these QPTs are in fact "crossovers" , 
because the two different kinds of physical behavior, su- 
pcrflow and soliton, can be connected by analytic con- 
tinuation. However, these QPTs have no meaning in the 
thermodynamic limit, where there is no Bose-Einstein 



condensation in ID, and as such are fundamentally re- 
stricted to the finite-size isolated systems typically found 
in experiments on Bose-Einstein condensates. Moreover, 
such crossovers can appear quite sharp in experiments, 
so that the matter of terminology becomes a question 
of theory, not experiment. In metastable states of mat- 
ter waves, such as soliton trains [l?! dl, the effects of 
dissipation can be suppressed and the metastable con- 
densate is observable. However, this picture does not 
extend into the medium- to strongly interacting regime, 
where quantum fluctuations cause mean-field solitons to 
decay [iSES]. Two questions follow. (1) Does the QPT 
indicated by mean- field analysis hold for stronger interac- 
tions? (2) If so, in what way is the system characterized 
on either side of this QPT, given that mean-field solitons 
are clearly no longer eigenstates? 

Our answer lies in the special class of many-body eigen- 
solutions called yrast states [U [22[ , defined as the lowest- 
energy solutions for fixed angular momentum. Studies of 
one-dimensional systems relevant to our chosen model 
have a long history, including exactly solvable quantum 
systems [a [H, HH, decay of persistent current (25H29| . 
and classical solitons [30l - [32l ]. In the thermodynamic 
limit it has been known since Lieb that there exist two 
excitation branches in the system, called type I and type 
II excitations. While the physical meaning of type I was 
clarified as the particle excitation and was found to agree 
with the Bogoliubov-type excitation in the weakly inter- 
acting limit, the meanin g of type II was elusive, described 
only as hole excitations [2j]. Seventeen years after their 
discovery, type II hole excitations were identified as a 
soliton branch by analysis of the energy of a classical soli- 
ton in terms of the nonlinear Schrodinger equation j30| . 
However, the validity range of the nonlinear Schrodinger 
equation is limited only within the range where the mat- 
ter wave possesses off-diagonal long-range coherence. 

The central finding of this article is that quantum 
solitons in the Lieb-Liniger Hamiltonian are precisely 
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the yrast states, and such states are the key to the 
metastable QPT previously identified in the mean-field 
context [lU [l(| . We first show how to distill the mean- 
field branches and QPT, which are previously found 
in the mean-field theory, from the metastable yrast 
states. Throughout the manuscript, this is our basic 
fashion of discussing the metastable states. The mean- 
field superflow-soliton QPT found in Rcfs. [HI, HH is 
shown to be obtained by extrcmizing the yrast spec- 
tra. In the weakly interacting regime, this type II ex- 
citation can indeed be called a soliton branch as shown 
in Ref. [30| because of quantitative agreement with the 
Gross-Pitaevskii mean- field theory [33[. We next intro- 
duce the concept of the "particle" and "hole" excita- 
tions which are clearly defined in the strongly interacting 
Tonks-Girardcau (TG) limit. We recover Lieb's result, 
and the metastable condition for the type II excitation 
branch leads to the observable quantum phase with a 
nonintegral single-particle average momentum. The type 
II "hole" excitation branch is made metastable (as op- 
posed to unstable) by subjecting the gas to a rotating 
drive, and is observable in typical "rotating-buckct-type" 
experiments [34j in a manner similar to the method used 
to create quantized vortices. Finally wc apply this con- 
cept to the regime of medium interaction strength. 

This article is structured as follows. In Scc.[II]we intro- 
duce the Licb-Linigcr Hamiltonian subject to an external 
rotating drive. The yrast problem, and basic properties 
of the eigenstates are described. In Sec. IIIII we investi- 
gate the many-body spectrum by exact diagonalization 
of the Hamiltonian in a truncated angular-momentum 
basis in the weakly interacting regime, comparing with 
those obtained by the mean-field theory. In Sec. IIVI wc 
study the opposite limit of the interaction strength, i.e., 
the strongly interacting TG limit where the many-body 
eigenproblem can be analytically solved using the Bosc- 
Fermi mapping. In Sec. |V] we address the intermediate 
regime of repulsive interaction between the weakly and 
strongly interacting limits via the finite-size Bethe ansatz 
approach. Finally, we summarize the results in Sec. I VII 



II. FORMULATION OF THE PROBLEM 
A. The model and yrast states 

We consider the same model as in Rcfs. [UGl, 

and 

solve its eigenproblem beyond the mean-field and Bo- 
goliubov theories. The Hamiltonian for periodic one- 
dimensional bosons with a contact interaction, 

N d 2 

3=1 3 3<k 

is known as the Lieb-Liniger Hamiltonian (LLH) Q, 
where 8j is the azimuthal angle that satisfies < 9 3 < 2ir, 
N the number of bosonic atoms, and <?id the effective 



strength of s-wave interatomic interaction in one dimen- 
sion (ID) [3i|. The length and energy units are the cir- 
cumference of the ring R, and ti 2 / (2mR 2 ) with m being 
the atomic mass, respectively. The coupling constant is 
measured in units of h 2 /(2mR) and c/id is hence dimen- 
sionless. The purpose of this article is to elucidate the 
many-body properties of these bosons subjected to a ro- 
tating drive. The LLH in a rotating frame of reference 
with an angular frequency 2f2 is given by 

H (fl) = H a - 2flL + n 2 N , (2) 

where 

j=l 3 

is the angular-momentum operator. From the singlc- 
valucdness boundary condition of the many-body wave 
function [34[ , one can show that solving the eigenproblem 
in the rest frame ifo^o = £'(0)\f r o suffices in order to ob- 
tain solutions to the eigenproblem H(fl)^f = E(il)^ [l6j ]. 
The eigenvalue is simply given by E(£l) = E(0) — 2Q(L) + 
£l 2 N, which is periodic with respect to f2. 

Throughout this article our approach is based on yrast 
problems [U [22[. Yrast, a Swedish term originally used 
in nuclear physics which can be translated as "dizziest," 
refers to the lowest energy state for a given angular mo- 
mentum. This approach is particularly profitable for a 
finite system, because all the information about physi- 
cal properties in the rotating frame are embedded within 
the spectrum in the rest frame. Thus the physical mean- 
ings of yrast states can be extracted by the simple trans- 
formation of yrast spectra. Since the LLH commutes 
with the angular-momentum operator, [Hq 7 L] = 0, the 
yrast problem is well defined irrespective of the sign and 
strength of interaction, and all the yrast states are eigen- 
states of both the Hamiltonians Hq and H(fl). 

All the eigensolutions of Hamiltonians Hq and H are 
classified according to the number of atoms N and total 
angular momentum L. Let us write the set of eigenstates 
classified into the subspace given by parameters (N, L) as 
\N, L; q) where q € |Z| is an energy quantum number that 
arranges the eigenvalues for fixed (N, L) in ascending or- 
der. The yrast states (the lowest-energy state under a 
given set of N and L) are denoted as \N,L;q = 1). The 
essential properties of the ground and low-lying excited 
states can be described within the yrast states. Thus 
we henceforth omit the quantum number q from the no- 
tations for eigensolutions. There are two external pa- 
rameters, the coupling constant and the external angular 
frequency of the rotating drive (divided by 2), written 
as (<7id,^), respectively. With the abbreviation of the 
quantum number q = 1 and for fixed coupling constant 
<7id, the eigenvalues that correspond to the yrast states 
are written as EVj^fi), where we explicitly write the 
parameter Q in the notation in order to clarify in which 
frame the system is. With this notation, the eigenso- 



3 



lutions in the rest (non-rotating) frame are written as 
E N , L (0). 

B. Center-of-mass rotation states 

Due to the translational invariance of the LLH with 
respect to 8 and f2, properties of a particular set of yrast 
states can be analyzed without solving the problem. We 
denote the set of yrast states for which total angular mo- 
mentum is equal to an integral multiple of the total num- 
ber of atoms N, as center-of-mass rotation (CMR) states. 
The energy of the CMR state takes the form 

E n ,l=jn(0) = J 2 N + V int , (4) 

where Vint is the interaction energy and J € Z is an inte- 
ger. We call J the center-of-mass quantum number, be- 
cause it physically expresses the amount of uniform trans- 
lation of the center-of-mass momentum. In the Gross- 
Pitaevskii mean-field theory, J is conventionally called 
the phase winding number; out of the mean-field regime, 
such terminology becomes questionable if not meaning- 
less. In the rotating frame, the energy of the CMR state 
is given by 

E N , JN (Q) = (J-n) 2 N + V int , (5) 

where the change in energy associated with the frame 
change is involved only in the kinetic energy term, and 
the interaction energy is completely separated from the 
parameter Q. 

For repulsive interactions <7id > 0, the ground state in 
the absence of the rotating drive is the state with zero 
angular momentum, -E/v,l=o(0). The excitation energy of 
the CMR states with a finite angular momentum L = JN 
is thus given by 

E N , JN (0) - E Nfi (0) = J 2 N, (6) 

which is independent of the strength of interaction guy . 
This is natural because changing the total angular mo- 
mentum by the amount JN is just a frame change and 
the interaction is isotropic. The ground state in the pres- 
ence of the rotating drive is characterized by the CMR 
quantum number 

Jo = [n + ^\, (7) 

where [^J denotes an integer that does not exceed x. 

Because of the periodicity in the eigensolutions, an 
cigenstate \N,L) with the energy En^(£1) has a de- 
numcrably infinite number of counterparts \N, L + JN) 
and En.l+jn(Q), corresponding to arbitrary values of 
J G Z. Solving the yrast problem for a limited range of 
fixed angular-momentum states, e.g., —N/2 < L < N/2, 
therefore suffices to obtain all the eigensolutions. More- 
over, the spectra are degenerate for the same magnitude 
of angular momentum, En,l — En,-l in the absence of 



rotating drive, while this degeneracy is resolved in the 
presence of rotation due to the Sagnac effect [3(|. All 
other yrast states for L out of this limited range can be 
obtained by shifting the total angular momentum by N 
while keeping the internal structure of the eigenstates. 
This is similar to a band theory concept, as discussed in 
Ref. [H], with -N/2 < L < N/2 playing the role of the 
Brillouin zone. 



III. WEAKLY INTERACTING LIMIT 

In our previous studies [HI, [lj|, we investigated the 
weakly interacting limit of the Bose gas on a rotating 
ring. The Gross-Pitacvskii equation, which corresponds 
to the mean-field approximation for the Hamiltonian (|2|) , 
has two kinds of solutions, namely, uniform superflow 
and soliton train [37| . The mean- field energy diagram 
is characterized by the set of bifurcations of the soliton 
branch from the superflow branch. These bifurcations 
make a continuous topological crossover in the conden- 
sate wave function possible by changing Q. The motiva- 
tion of this section is to demonstrate how in general to 
distill the mean-field branches from a sea of many-body 
eigenvalues. We argue how the mean-field soliton branch, 
for which average angular momentum is not quantized, 
emerges from the yrast spectra. The meaning of spectra 
related to symmetry breaking associated with the exis- 
tence of soliton branch is also discussed. 



A. Solution of the yrast problem 

To rewrite the LLH in second-quantized form, the 
bosonic field operator is expanded in terms of a plane- 
wave basis with the single-particle angular momentum 
I, 

, +oo 

H0) = -= Sie™, (8) 
V 2tt . — 

where the pre- factor of 1/v2tt comes from the normal- 
ization of the plane wave, and hi and SI are annihilation 
and creation operators which obey the standard commu- 
tation relations for bosons. Equation ((HJ) manifestly sat- 
isfies the periodic boundary condition i/j(0) = ip{9 + 2tt). 
The Hamiltonian ([TJ in second quantized form is then 
given by 

H = 

1 — — OC 

-foo 

+91D tyPlbrJ>nfa+l,m+n ■ (9) 

k, 1 ,771,71 — — OO 

Note that all the angular momenta are measured in units 
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FIG. 1: Yrast energy eigenstates for (a) repulsive, and (b) 
attractive interaction for N — 40 bosons on the ring, obtained 
by diagonalization of Hamiltonian ([9]) with the cutoff angular 
momentum \l c \ — 2. The spectrum has kinks where L is 
an integral multiple of N, and is symmetric with respect to 
L = 0, i.e., E n ,-l{Q) = En,l(0). In the large N limit, the 
number of points increases and the discrete yrast energies 
approach a continuous curve while the curvature and kink 
points remain unchanged. 

of h. The eigenstates can be expanded in terms of a Fock- 
state basis | {ni}) that represents the occupation number 
of each single-particle angular-momentum state, 

\{ril}) = | .. .,71-1,710,711, . . .) . (10) 

These states satisfy the conservation laws 

J2n t =N, Y, ln i= L - ( U ) 

( I 

In practice, for numerical calculations we require a cut- 
off angular momentum l c > 0. The range of the pos- 
sible total angular momenta for numerical diagonaliza- 
tion is hence limited to the interval of integer values 
L £ [—l c N,l c N]. In the weakly repulsive interacting 
regime giriN < 0(1), which we study in this section, 
a cutoff of l c = 2 pro vides a quantitative agreement in 
energy eigenvalues |38[ with those obtained by the Bethe 
ansatz shown in Sec. [Vj Thus all the results from this 
section are obtained with a cutoff of l c = 2. 

Figure Q] shows the yrast energies En,l{Q = 0) = 
(N, L\Hq\N, L), namely, the smallest eigenvalue obtained 
by the diagonalization of the Hamiltonian Hq within the 
restricted Hilbert space \N,L) for (a) <7id = 2.5 x 10 _2 7r 
and (b) g m = -2.0 x IO^tt for N = 40. The ratio 
of the mean-field interaction energy to the kinetic en- 
ergy corresponding to these values of <?id and N is (a) 
giT>N/(2ir) = 0.5, and (b) —0.4, respectively. The case 
of attractive interaction is shown just for reference, as 
our main interest is in the case of repulsive interactions. 
As the cutoff angular momentum is Z c ~ 2, AN + 1 yrast 
states (eigenstates corresponding to the eigenvalues for 
L £ {0, ±1, • • • ± 2iV}) are plotted. Recall that energies 
depend only on the magnitude of the angular momentum: 
E N , L = E N ^ L for O = 0. 

The key feature of the spectrum is that there appears a 
prominent kink at every CMR state L = JN. As shown 
in Sec. Ill Bl the excitation energy of these states is given 
by J 2 N. We note, however, that other states L ^ JN as 
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FIG. 2: (Color online) Energy eigenvalues of the Hamilto- 
nian ([2j obtained by the Legendre transformations (|12[) of 
yrast eigenstates in Fig. [1] Each curve is distinguished by a 
different total angular momentum L. The lower panels are 
enlargements of the upper panels near one of the swallowtail 
regions. The dashed curves show a comparison to the Gross- 
Pitaevskii superflow (blue) and soliton (red) branches. 

well as the curvature of the spectrum are also important 
and determine the existence of another quantum phase, 
as we show next. 

B. Superflow and soliton components 

In order to obtain the eigenstates of the LLH in the ro- 
tating frame, we transform the yrast spectrum according 
to the Legendre transformation, 

E N , L (0) ->■ E N . L {n) = E N , L (0) - 2flL + Q 2 N. (12) 

Figure [2] plots energies Sjy,i(fi) as a function of f2, where 
the finite number of dots, each of which is characterized 
by the different angular momentum L in Fig. [1] become 
convex downward curves in Fig. [2l Each curve is thus 
characterized by a different total angular momentum and 
has a minimum at a certain value of O. The degeneracy 
En,l(0) = £ , at,_l(0) in the absence of a rotating drive is 
resolved for finite O due to the Sagnac effect [H], i.e., the 
energy difference naturally arises in the corotating, and 
counter-rotating states with the external rotating drive. 

For repulsive interactions [Fig. HJa)], the energy 
En.l=j n{Q) corresponds to the ground state where 
Jo is the ground-state CMR quantum number given by 
Eq. ([7]). The angular- momentum states with L = JN 
correspond to the CMR states, and the center of the 
parabola is located at f2 £ {Z} at which the CMR state 
becomes the ground state. On the other hand, for at- 
tractive interactions [Fig. [2] (b)] the CMR state is not 
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FIG. 3: Energy eigenvalues of the Hamiltonian ([2]) that sat- 
isfy the extremization condition cIEl,n(Q) /dL = for TV = 40 
and (a) </id = 0.0257T and (b) (/id = — 0.027T. Each thin 
curve is characterized by an integral average angular mo- 
menta, L/N = J 6 {0, ±1, . . . }, while each thick point by 
different nonintegral angular momentum L/N ^ J so that 
each thick curve smoothly connects two thin curves with dif- 
ferent angular-momentum states. The energy of the former 
group agrees with the mean-field energy of the uniform su- 
perflow state, while the latter is well approximated by the 
mean-field energy of solitons. 

always a ground state, and is partially substituted for by 
the nonintegral average angular-momentum states. 

The transformation of the yrast spectrum according 
to Eq. (|12[) tells us that the eigensolutions of the Hamil- 
tonian in the rotating frame have an extremely high 
density of states around f2 € {±0.5, ±1.5, . . . } due to 
the crossing of many eigenvalues. These regions are en- 
larged in the lower panels of Fig. [2] for both signs of 
5id, where we also plot the energies of the stationary 
states given by the mean-field theory. Swallowtails were 
previously found within the mean-field theory to occur 
past the phase transition boundary between the uniform 
superffow and broken-symmetry soliton states fill ]. In 
the microscopic quantum theory, the high-density region 
also forms an upward/downward swallowtail-shape do- 
main for repulsive/attractive interactions, and the region 
is almost filled by various energy eigenvalues of various 
angular-momentum states crossing each other. The do- 
main with the high-density swallow-tail shape looks as 
if it is enclosed by the two kinds of stationary branches 
predicted by the mean-field theory. 

Although all the angular-momentum states shown in 
Fig. [2] are eigenvalues of the Hamiltonian ©, not all 
states are realized in practice. One example is vortex 
formation in a scalar condensate under rotation. Solv- 
ing the yrast problem in two dimensions results in all 
the angular-momentum states, including the rest con- 
densate (L = 0), off-axis vortex (0 < L < TV), a cen- 
tered vortex (L = TV), and vortex lattices (L > TV). In 
experiment, however, one drives the system with a spe- 
cific angular frequency. In such a situation, there exists a 
small distortion in the trap, which "selects" a metastable 
angular-momentum state with respect to the variation in 
the angular momentum of the condensate. As a result, 
in reality one does not observe a stationary off-centered 
vortex except as a transient state. 



(a) 




FIG. 4: (Color online) Average angular momentum L/N that 
gives dEL,N{Q)/dL = for (a) repulsive and (b) attractive 
interactions with TV = 40. Solid curves plot the single-particle 
angular momentum obtained by the mean-field theory. 

The same argument applies to our case. In the pres- 
ence of any kind of noise, such as an infinitesimal distor- 
tion of the trapping potential, quantum measurement of 
the matter wave, or whatever else breaks the translation 
symmetry of the ring trap, the realizable stationary state 
or metastable stationary state is determined by extrem- 
ization with respect to variations in angular momentum. 
In order to find the metastable states we impose the con- 
dition 

2%M = (13) 
oL 

with fl and <?id being fixed. 

Figure [3] plots energy eigenvalues that satisfy the con- 
dition (fT3j) as a function of tt; and Fig. |4] shows the 
corresponding angular momentum. These figures are 
quite similar to those given by mean- field theory, i.e., 
by imposing the stationary condition (|13[) for the mani- 
fold of eigenvalues we identify the mean-field stationary 
branches. The resultant branches are classified into two 
kinds according to the value of the angular momentum 
and have physical meanings as follows: 

Superflow: Due to the kink in the yrast spectrum at 
L = JN in Fig. [TJ these CMR states always satisfy the 
condition (|13[) . In particular, for the weakly interacting 
regime the CMR states can be specifically called uniform 



6 



supcrflow states, of which energies are given by Eq. J5]) 
and which correspond to the thin parabolic curves in 
Fig. El The energy of supcrflow states E^,jn agrees very 
well with the plane- wave energy 7V[(Q— J) 2 + gmN / (2tt)] 
in the mean- field theory pj| . 

Soliton Components: Other kinds of metastablc 
angular-momentum states appear that connect distinct 
supcrflow states as a function of fi, as shown by the thick 
curves in Fig. [3] The corresponding angular momentum 
divided by iV is nonintegral (see Fig.|4j, but it approaches 
integral values at both ends of this branch. These 
branches are equivalent to the maximum/ minimum en- 
velope of the high-density swallow-tail domain for repul- 
sive/attractive interactions, and can be approximated by 
soliton energies given by the Gross-Pitaevskii mean-field 
theory. 

However, we should not call the thick curve a soliton 
branch in a rigorous sense, because each point of this 
branch in Fig. [3] is the eigenvalue of the Hamiltonian and 
thus still possesses translational symmetry, unlike soli- 
tons. Instead we should call all the angular-momentum 
states inside the swallow tail in Fig. [2] the soliton com- 
ponents, because in the presence of infinitesimal noise 
these states do form a broken-symmetry state, which we 
denote |%). Soliton solutions of a Gross-Pitaevskii equa- 
tion can be interpreted in terms of the eigensolutions of 
the many-body Hamiltonian as a state where the several 
eigenvalues in the swallow-tail region are collectively su- 
perimposed. The derealization of a mean-ficld-like soli- 
ton in weakly interacting theories, as demonstrated by 
Dziarmaga et al. [39 |. is a dynamical demonstration of 
this idea. 

The energy associated with this superposition does not 
change significantly because the energy required to make 
it is on the order of 1/N. As a result, the energy of the 
broken-symmetry soliton state \x) is also well approxi- 
mated by the thick curve in Fig. [3] In the presence of an 
infinitesimal symmetry-breaking potential, the angular 
momentum is no longer a good quantum number. How- 
ever, the expectation value of the angular momentum 
(x\L\x) agrees well with that of the solitons obtained by 
mean-field theory, and thus behaves like that shown in 
Fig. H] [15]. With all these caveats in mind, we briefly 
say the branch drawn by thick curve in [4] is the quantum 
soliton branch in the weakly interacting regime. 

We also calculate the second derivative 
c1 2 En,l{&) / dL 2 with respect to Q, in order to check 
whether the metastablc angular-momentum state is a 
local maximum or minimum. For repulsive interac- 
tions the supcrflow state with a CM quantum number 
Jo = [O + 1/2J is indeed the ground state because 
the second derivative is positive at that point, while 
the thick points are local maxima with respect to L, 
since the second derivative is negative. For attractive 
interactions, the ground state is either a plane wave or a 
bright soliton in —0.5 < gmN / (2ir) < but the soliton 
becomes the sole ground state for gi£,N/(2ir) < —0.5. 
Consistently, the thick curve in Fig. [3] (b) becomes the 



global minimum. 

To sum up this section, we obtained the yrast states 
\N, L) of the LLH by diagonalization of the Hamiltonian 
in the weakly interacting regime. Among these eigenval- 
ues of the Hamiltonian in the rotating frame, we distilled 
the metastablc branches from the variety of yrast spec- 
tra by imposing an extremization condition. Two kinds 
of metastablc branches, superflow and quantum soliton, 
were found, consistent with mean-field theory. The re- 
gion where the different angular-momentum states in the 
quantum theory densely cross indeed agrees with the soli- 
ton regime predicted by the mean-field and Bogoliubov 
theories. The phrase "quantum flesh sewn onto classi- 
cal bones" has been used elsewhere [4l[ as a visual 
metaphor, perhaps inspired by x-ray images, to describe 
this accord. As we show later, the simple method shown 
in this section for obtaining metastablc states is appli- 
cable to other regimes; these two branches continuously 
exist over a wide range of interaction, from the weakly 
interacting regime all the way to the strongly interacting 
TG gas. 

IV. TONKS-GIRARDEAU LIMIT 

In the previous section, wc studied the weakly interact- 
ing limit of the LLH in the rotating frame to demonstrate 
how to obtain the mean-field-like stationary states. In 
the opposite limit of the strongly interacting TG regime, 
where the bosons are impenetrable and hence behave like 
spinless fermions, the eigenproblem can be calculated via 
the Bose-Fermi mapping. In this section we thus solve 
the yrast eigenproblem of free spinless fermions. In par- 
ticular, we introduce the particle and hole excitations, 
which are well defined in the fermionized gas, and show 
that these excitations are related to the mean-field sta- 
tionary states in the opposite weakly interacting limit. 

A. Bose-Fermi mapping 

The Bose-Fermi mapping theorem [l3| states that the 
eigenvalues Eb of impenetrable bosons are identical to 
those of spinless free fermions Ep of the same form of 
Hamiltonian, and the eigenfunctions of bosons $ b are 
generally written in terms of those of free fermions p 
as 

<S>B({e}) = *F({0})l[sgn(6 k -9 l ), (14) 

k>l 

where Ilfc>z s § n (^fc ~~ @i) is a unit antisymmetric function 
that takes the value +1 or —1 depending on the order 
of coordinates. This theorem holds for all the eigenso- 
lutions, and hence significantly simplifies our eigenprob- 
lem. The detailed properties of the TG gas are reviewed 
in Ref. 0. 

We first calculate the ground- and excited-state ener- 
gies of free fermions without taking the thermodynamic 
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limit. For simplicity of notation we show the analytic ex- 
pression only for an odd total number of particles. For an 
even number of particles the periodic boundary condition 
must be taken as antisymmetric. 

The ground state of N (odd) free fermions is obtained 
by the occupation of the lowest angular-momentum 
states from I = — lp to I = If [see L = in Fig. [5], 
where lp = (N — l)/2 is the Fermi momentum. The 
ground-state energy is thus 



If 

£jv,£=o(f2 = 0) = 22 l 2 = ^N(N 2 -l 



(15) 



l=-h 



We note that the N dependence of Eq. (fTSj) is the same as 
that of the bound state for attractive interactions, E oc 
—gioN(N 2 — 1), except for the prefactors (43j . 



£r^0 (ground state) 



- ti /, 



( I ) 

-h l F 
L=\ I I !******'* ! i I 



L=2 



I I****** 1 1 * 1 I 

; \ 

I I I *» 



(II) 

-!r If 

I I *»»»*» . I I. 



I I I * **** 1 * * I I . 



-h If 
-H- 1 *******! i I 



I, U 
*******! | . 




L=3M 




FIG. 5: (Color online) Low-lying yrast ground and excited 
states of free fermions. (I) Particle excitations where the an- 
gular momentum of a particle increases while a hole is po- 
sitioned at If + J- (II) Hole excitations where a particle is 
placed at the lowest unoccupied state and the angular mo- 
mentum of a hole decreases. When L is an integral multiple 
of N, only the center-of-mass angular momentum is shifted 
from the ground-state configuration. The lower panel shows 
type I and type II excitation energies as a function of L/N 
for N — 11 free fermions. 



B. Particle and hole excitations 

We next consider the low-lying excitations. Lieb has 
shown p3| that excitation of the repulsively interacting 
Bose gas in the thermodynamic limit has two branches. 
The first branch is called type I and was shown to be in 
agreement with the Bogoliubov spectrum of plane waves 
in the weakly interacting regime. The second branch is 
called type II, and this was supposed to be absent in 
the Bogoliubov spectrum. Intuitively, the type I and II 
branches correspond to the particle and hole excitations, 
respectively. 

We reconsider these branches in the context of yrast 
states. For the excited state Sjv,l(0) with total angular 
momentum L (7^ 0), there exist two kinds of excitations, 
type I (particle excitations) and II (hole excitations), as 
originally named by Lieb. To obtain these excitations 
one uses the following procedure. 

(I) Remove a particle at the Fermi momentum lp and 
place it at the momentum lp+L. For free fermions, there 
is no energy-level reconstruction in an (N ± l)-particlc 
system associated with removal or addition of a particle. 
The energy of the type I excited state E$ L (0) is thus 
obtained as 



e n]l(0) = E N ,o(0) -lp + (h + Lf 
Relative to the ground state the energy is 



(16) 



A£$ L (0) = E£> L (0) - E N , (0) =L(N + L- 1). (17) 

There is no limitation on the single-particle angular mo- 
mentum for this excitation. Such excitations are doubly 
degenerate for il — 0, for I — > — I. 

This type of excitation has an infinite set for the differ- 
ent CMR states J. The particle excitation of the L = JN 
state is achieved by removing a particle at lp + J and re- 
placing it at I = lp+J+L—JN. The resulting excitation 
energy is given by 



- p0) 



AE {1) 



n,l(°) = J N 
L > JN . 



+ (l F + J + L- JN) 2 - (l F + J) 2 , 

(18) 



(II) Starting from the ground state, remove a particle 
(create a hole) at the momentum lp — L + 1 and place 
the particle at lp + 1, where < L < N. It is clear from 
Fig. that the hole with this kind of low- lying excitation 
energy be created only within the range — lp < I < lp. 
The energy of this excited state is given by 

E%%(0) = E N , L (0) -(l F -L + l) 2 + (lp + l) 2 , (19) 

and the excitation energy is thus 

AE$ L (0) = L(N-L + 1), 0<L<N. (20) 

At L = N, the particle configuration in the angular- 
momentum space is the same as in the ground state, 
provided that the center-of-mass angular momentum is 
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shifted. Starting from the state L — N , we can consider 
the same kind of hole excitation, where we now place 
a particle at the lowest-unoccupied angular momentum 
If + 2, and place a hole at If + 2 — L + N. In this way, the 
type II excitation is extended for JN < L < ( J + 1)N 
states. This is a hole excitation of angular momentum 
lp — L + JN+ J+l with a particle fixed at lp + J, starting 
from the yrast state L = JN. Since the excitation energy 
of the CMR state L = JN is AE NtL = JN (Q) = J 2 N, the 
type II hole excitation energy is 

A<[(0) = E$ L (0)-E N , (p) 
= J 2 N+{l F + J + l) 2 

-(l F -L + JN + J+l) 2 , (21) 

where JN <L< {J +1)N. 

Excitation energies of type I and II are plotted in Fig. [5] 
for N = 11 for L £ {0, 1, . . . , 3N}, i.e., up to J = 2. 
The figure shows the spectra as continuous curves ac- 
cording to Eqs. (|17|) and (|21[) ; in fact the finite-size sys- 
tem discretizes these curves. We note that the region 
< L/N < 0.5 was presented in Ref. [H. 

C. Metastable hole excitation under rotation 




FIG. 6: Realizable yrast eigenstates in the strongly interact- 
ing TG limit for N = 11. All the energies are plotted relative 
to the ground-state energy of En,l=o(0)- The thin curves plot 
the energy of the states with angular momentum L — JN and 
L = ( J + 1)JV, while the bold curve is the type II excitation 
energy for JN < L < (J + 1)JV. The angular momentum 
decreases from (J + 1)N to JN according to Eq. (|23p along 
the bold curve. 



Next we consider the excitations under rotation, i.e., 
rotate all the yrast spectra according to J5jvi(f2) = 
E n ,l(0) ~ 2flL + n 2 N. The energy of the type II ex- 
cited state E$P L (Q.) for JN < L < (J + l)iV measured 
relative to En^ = q(0) is given by 



£^(O) = A<i(0)-2flL 



n 2 N. 



(22) 



In a similar manner to the previous section, we look 
for the metastable angular-momentum states by impos- 
ing the extremization condition dE^ L (Cl)/dL = 0. By 
inspection CMR states L = JN are (either ground or 
excited) metastable states. This extremization condition 
gives another metastable angular momentum, 



L = (N + 1)[J + 



1 



(23) 



and the corresponding energy, 



n-[J + - 



-I 2 



7V(iV+ 1) 



(24) 



As a function of fl this is a parabolic curve the mini- 
mum of which is located at 17 = J + 1/2 with energy 
N(N + l)/4, as shown in Fig. [5] Let us next compare 
this curve with those of two CMR states EN y L = jiy(£l) 
and E n _l=(j+i)n(Q)- These CMR energies intersect 
at = J + 1/2 with the energy E N j N (J + 1/2) = 
^iv,(j+i)Jv(^+l/2) = N/4. The minimum of ([24]) is thus 
higher than the value Em,jn by A^ 2 /4 at O = J + 1/2. 



Another important point is the emergence of certain 
critical angular frequencies where the metastable type II 
branch disappears and merges into the CMR branch: 



(25) 



The stable angular momentum approaches L = (J + 1)N 
at f2~, and L = JN at 17+, respectively, and the cor- 
responding energy coincides with the energy of CMR 
states. 

This is reminiscent of the uniform supcrflow to dark 
soliton transition in the weakly interacting limit, where 
there exists a critical angular frequency at which the soli- 
ton branch bifurcates from the superflow branch. Thus 
a continuous crossover between these topologically dis- 
tinct states exists in the TG limit, as well as the weakly 
interacting limit. Naturally, we can associate the hole 
excitations in the TG limit with the soliton branch in 
the weakly interacting limit. 



V. FINITE-SIZE BETHE ANSATZ 

We studied metastable states drawn from the yrast 
spectrum in the limits of the weakly and strongly inter- 
acting TG regimes by the extremization condition of the 
eigenvalues and by introducing the particle and hole exci- 
tations. The goal of this section is (i) to properly interpo- 
late these limits via the finite-size Bethe ansatz approach 
with the hint that the physical properties of the LLH with 
repulsive interaction are continuous and (ii) to show that 
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two distinct phases exist over the entire range of repul- 
sive interactions. We also vindicate our hypothesis that 
the soliton branch in the weakly interacting regime are 
"hole" excitations of the quasimomenta in general. The 
ground- and excited-state energies in the thermodynamic 
limit were given by the integral equation as a continuous 
limit of the Bethe equations [a, [23f . Recently, the spec- 
trum of LLH was obtained [44| by treating the inverse of 
the TG parameter, which is infinite at the TG regime, as 
the expansion parameter, and its analytical interpolation 
was given [45l ]. 



A. Bethe equations for quasi-momenta 

We first present how to derive the eigensolutions of 
the LLH by the Bethe ansatz 0. The delta- function 
interaction imposes the conditions 



d_ 



d 



_d_ 



d 



- 9mMe 1= e k ■ (26) 

Equation (f2l)|) is rewritten by the interchange of the sub- 
scripts j, k as 

d d \ J 1 



:.9id* 



8 i= 6 



(27) 



where ^ = ^({0}) is the many-body wave func- 
tion. The periodic boundary condition is expressed as 
*(0, 2 , • • • , On) = *(02, • • • , On, 27r). The entire coordi- 
nate space is expressed as 1Z : < 8-ji 1 < 0-r 2 < ••• < 
Otz n < 27r where {7?.i,7?.2, • • • ,T^n} is given by a permu- 
tation of {1,2,..., N}. The next procedure for solving 
the problem is to restrict the original coordinate space 1Z 
to the ordered space IZi and solve the Hamiltonian within 
the ordered space, say, IZi : < 6\ < ■ ■ ■ < 6^ < 27r. 
The LLH and the condition Eq. J27]) yield 



iY 



and 



00 



_d_ 

W 3 



^do 2 

3 = 1 3 



:5id* 



(28) 



(29) 



respectively, where refers to the many-body wave 
function in the ordered coordinate space. The periodic 
boundary condition and its derivative in the region IZi 
are given by 



§(0,1 



* n ) 



^(2tt,0 2 ,...,0n) 
V(6 2 ,..., N , 2tt), 



oe y 



7jV 



8=0 



'N : 



(30) 



(31) 



The Schrodinger equation ([28]) describes free particles. 
All eigenstates and spectra can therefore be represented 
formally by those of free particles. The eigenfunction in 
the region 1Z can be written as 



JV 



* = ]TA(7>,ft)exp i^vM) 



(32) 



where {£ n } are called quasi-momenta (or quasi-angular 
momenta). This is the basic idea of the Bethe ansatz: 
one writes down the many-body wave function in terms of 
a symmetrized superposition of plane waves with quasi- 
momenta, which implicitly includes all the effects of in- 
teractions. This wave function is a superposition of plane 
waves with N distinct quasi-momenta l\ < l 2 < ■ ■ • < 
£n, V means AH permutations of quasi- momentum in- 
dices, and A(P,1Z) is the coefficient of superposition 
of plane waves with a different configuration of quasi- 
momenta {£ n }. 

Substituting the wave function (|32j) into the conditions 
(|2T)1) and (f5D|) . we obtain the equations that determine the 
values of quasi-momenta {£ n }' 



where 



N -i2irli 



exp 



N 
L k=l 



Qkj = —2 arctan 



9lD 



j e {!,...,#}, (33) 



(34) 



is the two-body phase shift. 

Note that the quasi-angular momenta {£„} do not have 
a physical meaning per se; however, the sum of quasi- 
angular momenta does have a physical meaning as the 
total angular momentum, 



n=l 

Energy is also given in terms of {£ n } as 

E N , L (n = o) = ^2^, 



(35) 



(36) 



where the units of quasi-angular momentum and energy 
are the same as the previous sections, h and h 2 / (2mR 2 ), 
respectively. From the boundary condition we obtain N 
simultaneous nonlinear equations (Bethe equations), 



(-1) 



N-l^-2-nil 



JV 

n 



In - lm + ^fflD/2 
in - (-m - igiD/2' 



(37) 



which determine the set of values {£ n } for each atom 
n £ {1, . . . , jV}. Since all the quasi-momenta are known 
to be real and continuous for positive gm, Eq. (|37[) can be 
separated into real and imaginary parts, both of which 
are found to give the same set of solutions. Therefore 
it is sufficient to solve only the real part of the set of 
equations. 
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FIG. 7: (Color online) Comparison of exact yrast eigenener- 
gies of Eq. (|36|) obtained from the Bethe ansatz, yrast states 
obtained from the diagonalization of Hamiltonian ([9]) for 
N = 10, and the plane-wave energy branches for J £ {0, 1, 2} 
given by the mean-field theory. For Bethe ansatz and diag- 
onalization results, the ground (L — 0), and 2N yrast states 
\L\ 6 {1,. . . , 2iV} are plotted on a log scale. The inset en- 
larges the region 1 < <?id < 5 on a linear scale. 



B. Weakly interacting regime 

We numerically solve the real part of the Bethe equa- 
tions ([37]) for each set of energy levels characterized by 
the different total angular momenta. The numerical so- 
lution of Eqs. (|3T|) is highly sensitive to the initial set of 
trial values of {£ n }. If this initial set is sufficiently close to 
a solution for a target angular-momentum state, the set 
of solutions {£ n } can be correctly obtained. In contrast, 
if the initial set is closer to another angular-momentum 
state, the total angular momentum given by Eq. ([35]) re- 
veals undesired jumps, deviating from the target angular 
momentum. In such a case we again start from another 
initial set of trial values of quasi-momcnta. 

For simplicity we start with consideration of the trivial 
nonintcracting limit g^o = where all the quasi-angular 
momenta {£ n }, and hence the total angular momentum 
L, are zero for the ground state. The energy of the first 
excited state corresponds to the degenerate yrast levels 
-E/v,l=±i(^ = 0), where only one of the quasi-angular 
momenta has the value 1 or — 1 and the remaining quasi- 
angular momenta are zero. The second excited state has 
total angular momentum L = ±2, that is, two of the 
quasi-angular momenta take the value \£ n \ = 1. Higher 
excited states can be obtained in a similar way. Start- 
ing from these initial conditions, solutions can be ob- 
tained from the non-interacting to the strongly interact- 
ing regime by gradually increasing the value of giu . 

Results of the Bethe ansatz are obtained by the follow- 
ing steps: 

1. For L G {0, 1, . . . , [N/2\ }, the quasi-momenta 
are given by {£ n } G {(0,...,0), (1,0,..., 0), 



(1,1,0,. ..,0), (1,1,1,1,1,0,0,0,0,0)} for 

9id = and for N = 10, for example. Each 
set is directly calculated by solving Eqs. (f3"T|) with 
the target angular-momentum state from L = to 
L = N/2, respectively. 

2. The L G { LA72J , [N/2\ + 1, . . . , N} states are ob- 
tained from a transformation of the first [N/2\ — 
1 states of (i) via L = N — L, where L G 
{ [N/2\ , [N/2\ - 1,...,0}. The quasi-momenta 
are given by a transformation of the following 
form. For example, L = 1 with {£ n } = 
(1,0,0,0,0,0,0,0,0,0) transforms to L = 9 with 
{£ n } = (0,1, 1,1, 1,1, 1,1, 1,1). 

3. There exist degenerate spectra with the same mag- 
nitude of angular momentum El = E^l for the 
yrast states. The corresponding quasi-momenta are 
given by a transformation, L : {£\,l<i., . . . ,£isr) — > 
—L : (— £i, — £ 2 , ■ ■ ■ , —£n)- Although there exist a 
denumerably infinite number of other kinds of exci- 
tations of higher energy, these states are irrelevant 
for comparison with the previous results. Hence we 
do not use the larger quantum-number subscript 
q > 2, as in previous sections. 

4. We gradually increase giD from zero with the step 
size O(10- 4 ) for g 1D < 0(1), O(10- 2 ) for 0(1) < 
Aid < O(10), and O(K)- 1 ) for O(10) < g m < 
O(10 2 ). The convergence of the numerical solu- 
tions of Eqs. (|37[) is confirmed by comparing both 
sides of the Bethe equation with the substitution 
of the solution {£ n }. We set the tolerance factor, 
i.e., the difference in the left- and right-hand sides 
of Eqs. ([3T| . to be 10~ 8 . We also required as a sec- 
ondary convergence criterion that the total angular 
momentum be conserved to better than 10 -8 in the 
target angular momentum. If, during the changing 
of <?id, the angular momentum unexpectedly de- 
viates from the target angular momentum, and/or 
some of quasi-angular momenta show a jump as a 
function of gir>, these errors are detectable. For 
the attractive case, complex solutions of Eqs. (|3"T)) 
appear pfjj , indicating ground-state soliton forma- 
tion, which we do not treat here. 

In Fig. we compare low- lying excited states obtained 
by three different theoretical methods: Bethe ansatz, 
diagonalization, and GP mean-field theory. We note 
that the concept of yrast state for the angular momenta 
JN < L < (J + 1)N does not exist in the mean-field 
theory: this theory is concerned only with the single- 
particle angular momentum, which coincides with the 
average angular momentum in this theory. We thus plot 
the mean-field energy for the integral single-particle an- 
gular momenta. We plot the first AN + 1 yrast spectra 
En,l, \L\ G {0, . . . , 2iV}, as a function of the strength of 
interaction giD for ./V = 10. As expected, the rigorous 
Bethe ansatz spectra have the lowest energy for any giD , 
the mean-field plane-wave branch has the highest value, 
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FIG. 8: (a) Eigenvalues En,l of the Lieb-Liniger Hamilto- 
nian |l} obtained by the finite-size Bethe ansatz as a func- 
tion of interaction for L = (ground state), and L £ 
{±1, . . . , ±2JV} (excited states) where El = E-l- The hori- 
zontal dashed line corresponds to the energy of free fermions 
(|15|l with the same number of atoms, (b) Excitation ener- 
gies (En,l — En,o)/N of yrast states as a function of angu- 
lar momentum \L\/N for fixed strengths of interaction. A, 
B, C, D, and E correspond to the vertical line in (a). At 
L = JN G {0, ±N, ±2JV, ...} the energy is that of a CMR 
state (relative to the unboosted ground state), therefore inde- 
pendent of gin , and is given by J 2 , even in the thermodynamic 
limit. 



and diagonalization results have values in between those 
obtained by the Bethe ansatz and mean-field theory. For 
.9id ^ 1, the Bethe ansatz and diagonalization results 
quantitatively agree very well, while the latter becomes 
larger than the exact spectra for (/id > 1 [38j |- 



C. Medium- to strongly interacting regime 

We now further investigate the yrast states via the 
Bethe ansatz, going beyond the weakly interacting mean- 
field regime. Figure [8| a) plots En,l(P- = 0) for \L\ G 
{0, ±1, . . . , ±2iV} over a wide range of repulsive interac- 
tions, where the horizontal dotted line shows the ground- 
state energy of free fermions given by Eq. (fT5|) . All 
the ground- and excited-state energies monotonically in- 
crease with respect to <?id > 0. 

Note, however, that the energy does not monotoni- 
cally increase with respect to the total angular momcn- 
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FIG. 9: (Color online) Eigenvalues of the Hamiltonian ([2]) 
that satisfy the extremization condition 8El,n / dL = for 
(a) weakly interacting mean-field regime, (b) and (c) medium- 
interaction regime, and (d) strongly interacting TG regime. 
The corresponding average angular momentum of the thin 
curves is given by integers L/N = J. The thick curves, whose 
average angular momentum is noninteger, come from the type 
II excitation branch. 



turn for a fixed strength of interaction. This is clearly 
shown in Fig.[5](b), where we show the excitation energy 
(En,l — Emx=o)/N of yrast states for several values of 
fixed interaction strengths (indicated as A, B, C, D, and 
E) as a function of L/N. In the noninteracting limit, the 
the yrast spectrum is linear with respect to \L\/N with 
nodes at \L\/N G Z. While in a weakly interacting limit 
(plot A), the spectrum still looks almost linear, as the 
interaction g^o increases, the kinks in the yrast spectra 
at L = JN G {0, ±1, ... } become more pronounced due 
to the large increase in the energy of the yrast state in 
between JN < L < (J + 1)N. For strong interactions 
(curve E), the system is in the TG regime, which can be 
confirmed by the fact that the excitation energy has the 
value between Eq. ((2TJ) for JV = 9 and N = 11. 

Finally, we observe numerically that the excitation en- 
ergy of the CMR state L = JN is independent of giD , and 
is given by Eq. ©, namely {E N . L - E Nfi )/N = (L/N) 2 . 
This follows from the nature of the CMR state L = JN, 
which is just a Galilean boost of the nonrotating state; 
under this transformation interactions are unchanged. 

In order to see how these points are transformed in 
the rotating frame, we again rotate the yrast spectrum 
according to 

N 

i=i 

= £jv,l(0) -2nL + Q 2 N. (38) 
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FIG. 10: (a) Two CMR states with L = 0,N (thin curves) 
and the type II branch (thick curve) that connects them. The 
energy is defined relative to the interaction energy Vint- (b) 
The existence range of the type II branch with angular mo- 
menta L £ [0,N]. At the left (right) boundary, the angular 
momentum is L — N (L = 0), and L = N/2 at Q — 0.5 (ver- 
tical dashed line), decreasing linearly with respect to Q in the 
shaded area, (c) Enlargement of (b) in the weakly interacting 
regime. The solid curve is the phase boundary given by the 
Bogoliubov theory, (d) Difference between the phase bound- 
aries given by the Bethe ansatz and the Bogoliubov theory. 
The Bogoliubov phase boundary overestimates the extent of 
the shaded area. 



The results are shown in Fig. [5] for various strengths of 
interaction. The thin curves are parabolas (17 — J) 2 + Vint 
for various values of center-of-mass quantum numbers J. 
The lowest possible energy of the CMR state is thus given 
by VI„t at 17 = J e {Z}. 

The thick curves plot other metastable angular- 
momentum states, the angular momentum of which is 
given by a nonintegral multiple of N. The weakly in- 
teracting mean- field regime is shown in Fig. IHta), where 
the type II branch that satisfies the metastable condi- 
tion just starts to appear. Thus these are the energies of 
the quantum solitons [sec also Fig.^a)]. As the interac- 
tion increases [Figs. HJb) and IDJc)] the domain with the 
swallow-tail shape enclosed by the two CMR branches, 



as well as the size of the type II branch, increases. In 
the TG limit [Hid)], the area of the swallow-tail region 
saturates the spectra. 

These behaviors arc quantitatively summarized in 
Fig. [Tola), which shows the energy En,l{&)/N of 
metastable states relative to the interaction energy Vint at 
each strength of interaction. The CMR branches drawn 
by the thin curves no longer have a giu -dependence be- 
cause of the subtraction of Vint, while the thick curve 
gradually increases the domain over which it extends 
as giD increases. For simplicity we plot only two 
CMR branches with angular momenta L = 0,N, and a 
metastable state associated with the type II branch that 
smoothly connects these two CMR states. 

As 17 increases, the thick curve appears to bifurcate 
from the CMR branch with angular momentum L = N 
at a certain 17(< 0.5), and at 17 = 0.5 the energy becomes 
minimum. As 17 increases further, this branch smoothly 
merges into the CMR branch with angular momentum 
L = and eventually disappears at a certain 17(> 0.5). 
We therefore find that the same kind of energy bifurca- 
tion which was found in the mean-field theory persists 
over the full range of repulsive interactions. 

Figure fTDTb) plots the existence range of the 
metastable-state type II excitation branch. The shaded 
area indicates the existence of such a branch. For higher 
angular-momentum uniform superflow states L = ( J + 
l)N and L = JN, the boundary can be obtained by the 
parallel displacement of this figure along the 17 axis. The 
angular momentum on the phase boundary 17(< 0.5) is 
given by L = N, and the angular momentum linearly 
decreases as 17 increases, just like in Fig. 2] (a). At 
O = 0.5 (vertical dashed line), the value of the angu- 
lar momentum is given by L = N/2 irrespective of the 
strength of interaction. At a certain value of 0(> 0.5) 
the angular momentum eventually goes to zero, causing 
the metastable hole excitation branch to disappear. This 
behavior corresponds to the fact that in the mean-field 
theory the type II branch bifurcates from the plane- wave 
regime, developing nodes, and it again merges into the 
plane- wave regime with the increase of f2 flU [l6j . The 
phase boundary approaches f2 = (J + 1/2) ± A^/2 in the 
strongly interacting regime. 

In the lower panel of Fig. [TU] the phase boundary is 
compared with the one obtained by Bogoliubov theory 
in the weakly interacting regime. Bogoliubov theory thus 
predicts the quantitatively correct phase boundary to the 
5% level in the weakly interacting limit giu < 5 (for 
iV = 10), but it significantly overestimates the phase 
boundary as the interaction increases. 

These results indicate that the continuous change in 
the topologically distinct quantum phases can be found 
at any strength of interaction. For larger strength of 
interaction, the existence range of the type II branch in- 
creases. In this existence range the angular momentum 
changes linearly in f2, and the rate of change thus de- 
creases for larger coupling constant. 
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VI. CONCLUSIONS 

We addressed the continuous topological change in 
the repulsively interacting ID Bose gas on a ring, pre- 
viously found in the Gross-Pitaevskii mean-field the- 
ory [16[. In the mean- field theory the Gross-Pitaevskii 
equation has two kinds of solutions: uniform supcrflow 
and the broken-symmetry soliton train as a function of 
interaction strength and rotation. In the weakly inter- 
acting regime, the energy diagram is characterized by 
the smooth bifurcation of a soliton branch from a plane- 
wave branch in the rotating frame, which is the key to 
the continuous change in the topology of the condensate 
wave function characterized as a self-induced phase slip. 

In this article we vindicated this picture starting from 
the many-body Hamiltonian without assuming the exis- 
tence of the condensate wave function and spontaneous 
symmetry breaking. We solved the yrast problem of the 
original Licb-Linigcr Hamiltonian by three methods: di- 
agonalization of the Hamiltonian in the weakly interact- 
ing regime; Bose-Fermi mapping in the strongly interact- 
ing TG regime; and the Bethe ansatz approach for all 
regimes of repulsive interaction strength. 

We then obtained the eigensolutions in the rotating 
frame through transforming the eigenvalues according to 
specific values of the angular frequency of the external 
rotating drive Q. The extrcmization condition is im- 
posed so that eigensolutions which are realizable in prac- 
tice are extracted from a very large number of possibili- 
ties. The realizable states, namely those metastable un- 
der symmetry-breaking perturbations, reveals that two 
kinds of eigensolutions are physically distinguished. One 
is the superflow state in which angular momentum is an 
integral multiple of the number of atoms. The other is 
a quantum soliton characterized by a set of soliton com- 
ponents, which are also the yrast states. In the weakly 
interacting regime, the energy and angular momentum 
obtained by exact diagonalization and the Bethe ansatz 
agree well with those predicted by the Gross-Pitaevskii 
equation. This fact bears out the above physical mean- 
ings of metastable states. In the opposite limit, the 



strongly interacting TG limit was studied by the Bose- 
Fermi mapping. We introduced the concept of particle 
and hole excitations, which are well defined not only in 
the fermionized system but also for the whole interaction 
range in terms of the quasimomenta. The solution was 
similarly transformed into the rotating frame to extract 
the metastable states. 

In between the weakly interacting mean-field and 
strongly interacting TG limits, we employed the Bethe 
ansatz approach. In order to compare with the diago- 
nalization results, the set of Bethe equations was solved 
without substituting the summation with an integral, 
and the lowest (2N+1) discrete excited states were found. 
By the same transformation into the rotating frame, we 
elucidated how the energy diagram of these topologically 
distinct states changes as the strength of interaction in- 
creases. 

Energy and angular momentum of the two kinds of 
topologically distinct states exist over the whole of repul- 
sive interaction guy . The quantum phase diagram in the 
f2-<7iD plane for the quantum soliton with a single den- 
sity notch and with angular momentum < L < N was 
explicitly shown. This metastable quantum phase transi- 
tion is technically a crossover: it can occur only in a finite 
system, as expressed by our choice of units in terms of 
the ring circumference, and all states are connected ana- 
lytically. Nevertheless, one finds a sharp change between 
distinct physical states which will appear as a QPT in 
experiments. 
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